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- - - Self-consistent system of spinor, scalar and BI gravitational fields in presence of 

■ magneto- fluid and A-term is considered. Assuming that the expansion of the BI 

. universe is proportional to the a/ component of the shear tensor, exact solutions 

for the metric functions, as well as for scalar and spinor fields are obtained. For 
a non-positive A the initially anisotropic space-time becomes isotropic one in the 
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X/^ • process of expansion, whereas, for A > an oscillatory mode of expansion of the BI 

^ . model occurs. 
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m ; 1. INTRODUCTION 

o 

The discovery of the cosmic microwave radiation has stimulated a growing interest 
I . in anisotropic, general-relativistic cosmological models of the universe. The choice of 
^1 anisotropic cosmological models in the system of Einstein field equation enable us to study 
the early day universe, which had an anisotropic phase that approaches an isotropic one 
Bianchi type I (BI) cosmological models that are anisotropic homogeneous universes play 
an important role in understanding essential features of the universe, such as formation of 
' galaxies during its early stage of evolution. An LRS BI model containing a magnetic field 
directed along one axis with a barotropic fluid was investigated by Thorne [2| . Jacobs 0, |^ 
investigated BI models with magnetic field satisfying a barotropic equation of state. Bali 
studied the behavior of the magnetic field in a BI universe for perfect fluid distribution. 
In this paper we study the self-consistent system of spinor, scalar and BI gravitational 
fields in presence of magneto-fiuid and cosmological constant. Solutions of Einstein equa- 
tions coupled to a spinor and a scalar fields in BI spaces have been extensively studied 
by Saha and Shikin 0, Q S 0. In the aforementioned papers, we considered spinor 
field in BI universe where nonlinearity occurred either due self-coupling or induced by an 
interacting massless scalar field. Here, considering the BI universe filled with magneto- 
fiuid we make an attempt to study a system, where all the four fields, i.e., scalar, spinor, 
electro-magnetic and gravitational ones, play active part in the evolution process. 
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2. FUNDAMENTAL EQUATIONS AND GENERAL SOLUTIONS 

We choose the action of the self-consistent system of spinor, scalar and gravitational 
fields in the form 

y{g;XI/,Xj/,(p)= I {R + ^)^da, (2.1) 

where R is the Ricci scalar and L is the spinor and scalar field Lagrangian density chosen 
in the form[^ 

-mi^i/A+l(p,«(p'«(l+AF). (2.2) 

Here A is the coupling constant and F is some arbitrary functions of invariants generated 
from the real bilinear forms of a spinor field. We choose F to be the function of / = = 
(YV)^ and J — P^ — (ixffY^ Xj/)'^ , i.e., F = F{I,J), that describes the nonlinearity in the most 
general of its form [8|. As one sees, for A = we have the system with minimal coupling. 
The gravitational field in our case is given by a Bianchi type I (BI) metric in the form 

ds^ = alidx^f - a\{dx^f - al{dx^f - ajidx^)^, (2.3) 

with ao = \, x^ = ct and c = 1. The metric functions Ui {i = 1,2,3) are the functions of 
time t only. 

Variation of 1)2.11) with respect to spinor field y/{^f) gives nonlinear spinor field equa- 
tions 

iY^y^\j/-m\j/+&\]/ + ^i)^\]/ = 0, (2.4a) 
iy^\[fY^+m\[r-^\{f-^i\{fy^ = 0, (2.4b) 

where we denote 

^ = XS(p^a(P'"dF/dI, ^ = XP(p^a(p'''dF/dJ, 

whereas, variation of ()2.1|) with respect to scalar field yields the following scalar field 
equation 




Varying 1)2.1)) with respect to metric tensor g^y one finds the gravitational field equation 
which in account of cosmological constant A has the form 
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Here K is the Einstein gravitational constant and over-dot means differentiation with 
respect to t. The energy-momentum tensor of the system is given by 

= ^g^''(^n^v¥+¥rv^fi¥-^^i¥rv¥-^vnv^ (2.7) 

The energy-momentum tensor of the magneto-fiuid is chosen to be 

T;im) = i^ + P>^^"-pS",+E;, (2.8) 
where E^y is the electro-magnetic field given by Lichnerowich 

E;^11 [\h\^{u^u' - ^5;) - h^h''] . (2.9) 
Here u^^ is the flow vector satisfying 

g^yu^u'' = h (2.10) 
jU is the magnetic permeability and hjj, is the magnetic flux vector defined by 

/^M = -^*i^vM"^ (2.11) 
where ★F^v is the dual electro-magnetic field tensor defined as 

•^F^iv = ^'^-^^livapP"'^ ■ (2.12) 

Here F^^ is the electro-magnetic field tensor and is the totally anti-symmetric Levi- 

Civita tensor with £0123 = +1. Here the comoving coordinates are taken to be = 1, = 
= u^ = 0. We choose the incident magnetic field to be in the direction of x-axis so that 
the magnetic flux vector has only one nontrivial component, namely h\ ^Q. \a view of 
the aforementioned assumption from (j2.1ip we obtain F12 = F13 = 0. We also assume that 
the conductivity of the fluid is infinite. This leads to Fqi = Fq2 = F03 = 0. Thus we have 
only one non- vanishing component of i^tv which is F23. Then from the first set of Maxwell 
equation 

Ffiv;P+fvl5;^i+P'l3^i-v=^^ (2.13) 
where the semicolon stands for covariant derivative, one finds 

F23 = ^, ^ = const. (2.14) 

Then from (j2.1H) in account of (j2.12j) one finds 

h, = ^. (2.15) 
IJ,a2a3 

Finally, for Ej^ we find the following non-trivial components 

£0 = El = -El = -El = (2.16) 
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In fl2.8|) £ and p are the energy and pressure of the fluid, respectively. In this note we 
assume that the equation of state 

p = i;e (2.17) 

holds. Here C, varies between the interval < ^ < 1, whereas <^ = describes the dust 
Universe, C = j presents radiation Universe, ^ < ^ < 1 ascribes hard Universe and 
^ = 1 corresponds to the stiff matter. The Dirac matrices y^{x) of curve space-time are 
connected with those of Mincowski space as 

f = f/a^, y^ = ya^, Ai = 0,1,2,3. (2.18) 

In the Eqs. ()2.4p and ()2.7p is the covariant derivatives acting on a spinor field as 



111 

^M^ = ^-rMV^, v^v^=0 + v/r^, (2.19) 



where are the Fock-Ivanenko spinor connection coefficients defined by 

T^ = \f{Tl^yv-d^ya). (2.20) 
For the metric ()2.3p one has the following components of the spinor connection coefficients 



T^, = {l/l)d^tf. (2.21) 

We study the space-independent solutions to the spinor and scalar field Eqs. ()2.4|) and 
(|2.5p so that i//= !//(?) and (p = (p{t). defining 

T = aoaia2a3 = (2.22) 

from ()2.5|1 for the scalar field we have 

(p = C j[x{\+XF)Y^dt. (2.23) 

Setting Vj{t) = ^\irj{t), j =1,2,3,4, in view of (jTTTI|l and ^^TI^ from (jTIal one 
deduces the following system of equations: 

Vi + i{m-&)Vi-mi = 0, (2.24a) 

V2 + iim-^)V2-':^V4 = 0, (2.24b) 

V3-i{m-^)V3 + mi = 0, (2.24c) 

V4-i{m-^)V4 + ':^V2 = 0. (2.24d) 

From ()2.4a|) we also write the equations for the bilinear spinor forms 5, P and A'^ = 

5o-2$fA[j = 0, (2.25a) 

/b-2(m-^)A|] = 0, (2.25b) 

A{] + 2(m-^)/'o + 25fSo = 0, (2.25c) 



where Qq = tQ, leading to the relation S-^+P^ + (A^)^ = C^/ T^, = const. As one sees, 

for F = F{I) ()2.25ap gives S = Cq/t, whereas for the massless spinor field with F = F{J) 
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()2.25b|l yields P = Dq/t. In view of it for F = F{I) we obtain the following expression for 
the components of spinor field 

(2.26) 

with Q being the integration constants and are related to Cq as Cq = + C2 — — C^- 
Here j8 = J{m — &)dt. In case of F = F[J) for the massless spinor field we get 

(2.27) 

1^3 = T-i/2(/Die'^+D3e-'^), v/4 = T^^/^ (^-^^e'a +^4^-'^) . 

The integration constants D,- are connected to Dq by Dq = 2 {D\ +D2 —D^—D^). Here we 
set o = J'^dt. 

Once the spinor functions are known explicitly, one can write the components of spinor 
current = xffy^y/, the charge density of spinor field p = {jo- j^)^^'^, the total charge of 

00 

spinor field Q= J p ^/—^gdxdydz, the components of spin tensor 5^^'^ = |v^{y^C7^^ + 

— 00 

a'^^y^ji// and other physical quantities. 

Let us now solve the Einstein equations. In doing so we first write the expressions for 
the components of the energy-momentum tensor explicitly: 



mS + C^/2T^{\+XF)+e + ^^^^, (2.28a) 



tI = &S + '^P-C'^/2%'^{\+XF)-p + (2.28b) 



Ti = ^S^^P-C^llx\\+XF)-p (2.28c) 



tI = ^5 + ^P-CV2t^(1+AF)-;? (2.28d) 

(2.28e) 

In view of = 7^3^ from ()2.6b|l . (j2.6c|) we find 

a2 = a3Dexp(x J (2.29) 

with the constants of integration D and X being integration constants. 

Following Bali let us assume that the expansion (0) in the model is proportional to 
the eigenvalue a/ of the shear tensor o^. Since for the BI space-time 
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the aforementioned condition leads to 

ai = (a2a3)^, (2.32) 

with A'^ being the proportionahty constant. 

In account of ()2.22j) from ()2.29p and ()2.32j) after some manipulation for the metric 
functions one finds 

ai = T^/(^+i), (2.33a) 

a2 = v/DTi/2(^+i)exp[| /-], (2.33b) 



as = ^Ti/2(^+i)exp 
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As one sees from ()2.33|) for z ^ t" with n > I the exponent tends to unity at large t. 
In this case the anisotropic model becomes isotropic one iff Z) = 1 and A'^ = 1/2. Let us 
also write the invariants of gravitational field. They are the Ricci scalar Ii = R ^ I/t^, 

h = R^vR^^ = RjiRy ~ 1 /t^"^ and the Kretschmann scalar I3 = Rap^yR"^^^ ~ 1 /t^"^- As we 
see, the space-time becomes singular at a point where T = 0, as well as the scalar and 
spinor fields. Thus we see, all the functions in question are expressed via T. In what 
follows, we write the equation for r and study it in details. 

Summation of Einstein Eqs. ()2.fia|l . ()2.6bjl . ()2.fjc|l and ()2.6djl multiplied by 3 gives 

f 3 / 2J^2 X 

-^-K{mS + ^S + WP + e-p + — To -3A. (2.34) 

For the right hand side of Eq. ()2.34j) to be a function of T only, the solution to this 
equation is well known In what follows we study this equation for some concrete 

form of F. In doing so let us demand the energy-momentum to be conserved, i.e., T^.y = 0, 
which in our case takes the form 

luTo^y _ ^r/ _ - ^T^ = 0. (2.35) 
T ^ " ^ ai a2 (33 ^ 

In account of the equation of state ()2.17|) and 

(m-^)5o-^/b = 

which follows from ()2.25|1 . after a little manipulation from ()2.35|1 we obtain 

e = eo/Ti+^, p=^eo/T'+^. (2.36) 

Let us recall that we consider F as a function of I, J or I ± J. If we choose F = F{I), 
then setting m = we come to the analogical equation corresponding to the massless 
spinor case with F =F{J) or F = F{I±J), whereas setting A = we have the system with 
minimal coupling. Under this assumption from ()2.25|) one finds 

5=^. (2.37) 
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In view of flTTTjl . (jT^ and the fact that a2a3 = T^/(^+^\ Eq. (1^1) can be rewritten as 

f = ^ K(mCo + ^Co + eo( 1 - + - 3A. (2.38) 

Recalling the definition of S) the solution to Eq. ()2.38|) can be written in quadrature 

J K(mCoT + C2/2(l + AF) + eoTi-^ + ((A^+ l)^2/3/iA^)T2^/(^+i)) - At2 + £ 

(2.39) 

with E being some integration constant. It should be mentioned that being the volume- 
scale T is non-negative. At the points where T = there occurs space-time singularity. On 
the other hand, the radical in ()2.39p should be positive. This fact leads to the conclusion 
that for A > the value of T is bound from above as well, giving rise to an oscillatory 
mode of expansion of the BI universe. For a non-positive A, we have picture with fast 
expanding x with time. 

3. CONCLUSIONS 

A self-consistent system of spinor, scalar and gravitation fields has been studied in 
presence of magneto-fluid and cosmological term A. With the presence of F23 component 
of electro-magnetic field tensor, the system can be viewed as one where all the four fields, 
i.e., scalar, electro-magnetic, spinor and gravitational ones, are taken into consideration. 
Assuming that the expansion of the BI space-time is proportional to the <3\ component of 
the shear tensor, solutions for the metric functions a,(?) are obtained explicitly in terms 
of volume-scale x. Expressions for the scalar and spinor fields are also obtained in terms 
of X. For the non-positive A we obtain exponentially expanding BI universe, which means 
the initially anisotropic space-time becomes isotropic one in the process of expansion. 
For a positive A an oscillatory mode of expansion takes place. Choosing the integration 
constant E and initial value value of x it is possible to obtain solutions those are regular 
everywhere. A detailed numerical study of the Eq. (|2.38|) we plan to perform in short. 
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